Small-q Phonon Mediated Unconventional Superconductivity in the Iron Pnictides 
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We report self-consistent calculations of the gap symmetry for the iron-based high-temperature 
superconductors using realistic small-q phonon mediated pairing potentials and four-band energy 
dispersions. When both electron and hole Fermi surface pockets are present, we obtain the nodeless 
s± state that was first encountered in a spin- fluctuations mechanism picture. Nodal s± as well as 
other gap structures such as d x 2_ y 2, s± + d x 2_ y 2 and even a p-wave triplet state, are accessible upon 
doping within our phononic mechanism. Our results resolve the conflict between phase sensitive 
experiments reporting a gap changing sign attributed previously only to a non-phononic mechanism 
and isotope effect measurements proving the involvement of phonons in the pairing. 
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One of the foremost issues in contemporary condensed- 
matter physics is the nature of the medium-high- 
temperature superconductivity (SC) in the recently dis- 
covered iron pnictides The focal point of current re- 
search is to understand whether the bosonic pairing that 
drives the SC is due to spin- fluctuations or phonons, and 
how the pairing mechanism concurs with the emerging 
symmetry of the superconducting gap [2]. Experimen- 
tal techniques have recently given controversial results 
0, Q • Angular-resolved photoemission spectroscopy and 
Andreev spectroscopy indicate one nearly isotropic gap or 
two isotropic gaps (EH?} , whereas penetration depth [H, [s| 
and nuclear magnetic resonance (NMR) To|, II] measure- 
ments support both nodeless and nodal gap structures. 
Moreover, very recent experiments on the evolution of the 
SC gap symmetry of LaFeAsOi-^F^ with doping, con- 
clude that there is a transition from nodeless to nodal 
SC upon moving from the under doped to the over doped 
regime of this compound [12[ . Notably, a nodeless gap is 
not necessarily a conventional isotropic s- wave gap here. 

An important ingredient for discussing the gap symme- 
try is the Fermi surface (FS) of both 1111 (LaOFeAs) and 
122 (BaFe2As2) parent compounds, which is remarkably 
simple. The strongly two-dimensional FS mainly con- 
sists of two electron and two hole pockets, with a nesting 
wave vector Q = (tt^tt) connecting them in the folded 
Brillouin zone (BZ) [13|, LL4[ Fermi surface nesting stabi- 
lizes a spin density wave which yields to superconduc- 
tivity upon doping jTjjj. If spin- fluctuations mediate the 
pairing then a nodeless unconventional superconducting 
gap reversing its sign between the electron and hole FS 
sheets is expected, known as the s± state @, Q3, EH, E3 • 

Phase sensitive experiments are needed to distinguish 
the s± state from the conventional isotropic s- wave state 
and indeed, inelastic neutron scattering results appear to 
be compatible with the coherence factors of the s± state 
This has been interpreted as evidence for a non- 
phononic pairing mechanism probably involving spin- 
fluctuations. However, a systematic NMR study of the 



dependence of low-lying excitations upon fluorine dop- 
ing in LaFeAsOi-^F^ has shown that antiferromagnetic 
(AFM) fluctuations are suppressed with doping while T c 
remains unaffected, imp lying that AFM fluctuations may 
be irrelevant to SC [18| . Most importantly, isotope effect 
measurements Q directly prove the involvement of the 
phonons in the pairing mechanism. Additional experi- 
mental findings do support the relevance of phonons in 
the pairing |l9[. However, it is mysterious how this could 
be compatible with a gap that changes sign. 

In this Letter we propose a new picture for super- 
conductivity in the iron pnictides on the basis of self- 
consistent calculations of the gap symmetry. We argue, 
first, that small-q electron-phonon scattering dominates 
the pairing. Second, we show that this solely phononic 
mechanism in the presence of strong Coulombic interac- 
tions can naturally produce the unconventional gap sym- 
metries proposed for the pnictides. When both electron 
and hole FS pockets are present, the nodeless s± gap is 
obtained self-consistently, because of the intricate four 
band pnictide FS combined with strong interband scat- 
tering. Third, nodal gap structures may also develop 
when going from underdoped to overdoped regimes upon 
doping. These include s±, d x 2_ y 2-wa,ve, s± +d x 2_ y 2 and 
even p-wave triplet SC. 

The effective electron-phonon interaction may be dom- 
inated by the forward processes only when the Coulomb 
interaction is short range or Hubbard-like 20|, 21], which 
is a requirement for antiferromagnetism as well. There- 
fore, the proximity of the AFM and small-q phonon me- 
diated SC phases is natural, and this is indeed the case in 
the pnictides. Unconventional SC due to small-q phonon 
pairing has already been considered for high-T c cuprates 
j2ll-[23|, heavy fermion 24], organic [IBj and cobaltite (26| 
systems. 

Our aim is to describe the electronic structure of the 
iron pnictides around the Fermi energy as realistic as pos- 
sible. To this end we construct a four band tight-binding 
(TB) model which captures the essential physics of the 
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FIG. 1. (Color online) a) Fermi surface contour of the 
undoped pnictide compounds in our four-band model, b) 
Wavevector dependence of the static spin susceptibility along 
high symmetry lines. The peak indicates well nested electron- 
hole pockets at Q = (tv,tv). 



undoped compounds, while combining features from the 
previously introduced TB models [27j . The energy bands 
are: 

e\ 2 = — 1\ 1,2 (cos k x + cos k y ) — t\ cos k x cos k y — e^ 1 ' 2 , 



-1,2 



—t\ (cos k x + cos ky) — t\ cos k x cos k y 
±£3 sin k x sin fc^ , 



with t\ 



-1.14, 



-0.2, t§ = 0.2, ef. 



(1) 
-1.5, 



q 1 = -0.3, q 2 = -0.2, t§ 



0.24, e^ 1 = 0.6 and 
e^ 2 = 0.4. We have verified that these TB dispersions 
are in good agreement with band structure calculations; 
not only the FS topology and the susceptibility are repro- 
duced (see Fig. [T]) but also the angularly resolved density 
of states (arDOS) at any Fermi momentum kp is properly 
given (arDOS ~ |Ve^2 (kp)| -1 ). To probe the doping 
dependence of the gap symmetry in our self-consistent 
calculations, we have introduced a chemical potential \i 
varying in the interval —0.6 < \i < 1.4. When \i > 
{fi < 0), the system is electron (hole) doped. The Fermi 
surfaces at various chemical potentials as well as the cor- 
responding arDOS are depicted in Fig. [2l 

We solve the most general multiband Bardeen-Cooper- 
SchriefTer (BCS) equation at zero temperature [28|, 



A„(k) = -^ VW(k,k') 



A ra .(k / ) 
2£„<(k') ' 



(2) 



where A n (k) is the momentum dependent superconduct- 
ing gap and E n {k) = ^e^(k) + A^(k) the quasiparticle 
excitation spectrum of the n th band. The pairing poten- 
tials y nn /(k, k') can support intra (n = n') and interband 
(n ^ n') scattering. As a first approximation, we will 
take VW(k,kO = V(k,k / ). This yields a "global" SC 
gap symmetry valid over the entire B.Z., which is nev- 
ertheless driven by the characteristics of each seperate 
band, such as the arDOS distribution. 

The effective pairing potential of the small-q phonon- 
mediated interaction takes the form [23] : V(k, k') = 
^c* — g 2 + |^ k /|2 properly continued periodically, where a 
repulsive Coulomb pseudopotential V* is responsible for 



screening at short distances whereas the negative term is 
the attractive phonon part. The pairing kernel is charac- 
terized by a smooth momentum cutoff q c which selects 
the small wave vectors in the attractive phonon part 
while at larger wave vectors the repulsive Coulomb pseu- 
dopotential may prevail. Decreasing the cutoff q c leads to 
a situation that has been named momentum decoupling 
(MD) [2l|. In the MD regime there is a tendency for su- 
perconductivity to decorrelate in the various FS regions; 
the gap function gradually loses its rigidity in momentum 
space. In regions of the FS with high arDOS we observe 
higher gap amplitudes, whereas in the low arDOS regions 
the gap is smaller. This arDOS driven anisotropy of the 
gap is the fingerprint of MD [2l| and may be behind mul- 
tiple gap signatures in the experiments. 

To solve the BCS equation self-consistently with the 
momentum dependent pairing kernel, we perform the 
convolution integral in Eq. ((2]) by utilizing an FFT cyclic 
convolution technique. The k-summation was performed 
over a512x512 grid in a thin shell around the FS, in order 
to get maximal resolution. In the results presented here, 
the Coulomb pseudopotential is fixed at V* = 0.09V p h 
and the pairing amplitude at V p h = 2. Several values of 
the cutoff parameter q c were considered varying from tt 
to 7r/8 and various values of the chemical potential were 
explored. 

When q c > 7r/5, anisotropic s- wave, sign-preserving 
solutions are found at all dopings. Further reducing q c 
yields interesting doping-induced transitions between dif- 
ferent gap symmetries. When 7r/5>g c >7r/8we get re- 
markable self-consistent solutions of the gap, depending 
on the doping level, that include nodal s±+d x 2_ y 2-wdive, 
s±, d x 2_ y 2-wave and p- wave symmetries. Here we present 
calculations for q c = 7r/6 as a representative example of 
our results. Fig. [2] depicts the evolution of the FS, colored 
by arDOS -1 as a function of doping. In the same plots, 
the lines where the respective self-consistently calculated 
SC gap vanishes are drawn. 

Starting from the electron doped side, we find an 
anisotropic d x 2 _ y 2 -wave gap in the regime 1.4 > \i > 1.0 
as shown in Fig. 0k . I n this doping region, the high-DOS 
hole pockets lie away from the Fermi level, and only the 
electron pockets drive superconductivity. The SC gap 
contains nodes intersecting the FS in the lowest arDOS 
points (Fig. [2^). The arDOS distribution at these dop- 
ings is highly anisotropic with the points where the gap 
nodes intersect the FS having considerably lowest DOS. 
Thus, gapping the rest of the FS while having nodes at 
the lowest DOS points is energetically favorable for SC. 

Decreasing \i shrinks the electron pockets and makes 
the arDOS more isotropically distributed. When 0.95 > 
ji > 0.65, a structure of the p-wave (sin k x ) type develops 
with nodes at the locations where the electron pockets 
cross (Figs. [2] and |Ud). Although these FS points have 
maximal arDOS this is only approximately two times 
larger than the minimum arDOS on the FS. The sys- 
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FIG. 2. (Color online) Contours of the FS, in the folded BZ, 
colored with arDOS -1 for u = (1.0, 0.75, 0.45, 0.1, -0.4, -0.6) 
from panel a) to panel f), respectively. Red FS lines means 
small arDOS. In the same panels we plot with blue lines the 
points where the respective, self-consistently obtained, SC gap 
vanishes (q c = ty/6). Around half filling we find an s± gap 
that persists in the electron (/x = 0.1 (d)) as well as hole 
doped region (see FIG. [4b). The s± gap is nodeless with its 
vanishing area almost halfway between the electron and hole 
pockets and persists at further electron doping even when 
the hole bands do not contribute to the FS (c). More electron 
doping results in a cascade of transitions to triplet (sin fee) (b) 
and (^,2 _ y 2 -wave (a) solutions. For \i < 0, the nodes of the 
s± solution move closer to the outer hole pocket and finally 
intersect it (e) . In the doping limit when the FS becomes like 
in (f) the SC gap is a superposition of s±+d x 2_ y 2 symmetries. 
Note that the hole FS is thicker than the electron one due to 
its higher DOS, therefore we have chosen to use a different 
color grade for the electron and hole bands. 



tern can still minimize its energy by keeping the rest of 
the FS gapped while leaving nodes on these points. Since 
superconductivity in our approach originates from parti- 
cle pairing in the same band, the symmetry of the gap 
function implies a spin triplet order parameter. Such a 
state has already been proposed for the iron pnictides in 
a different context [29[. However, a spin-singlet p-wave 
order parameter is also possible if one considers interband 
pairing (30| . 

In the interval 0.6 > /i > 0.25 a nodeless s± gap is 
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FIG. 3. (Color online) A typical, self-consistently obtained 
sign changing s± superconducting gap, originating from a 
small-q phonon interaction. The zeroes of this gap emerge 
between the electron and hole pockets (for q c — tt/6 and 
ji = —0.1). Notice the anisotropy of the gap function and 
how the gap amplitude varies over different regions 



found, yet with the nodal lines close to the electron pock- 
ets (Fig. [2b). The anisotropy of the arDOS is further re- 
duced and any node on the FS is not the preferred con- 
figuration. This s± solution is mainly a superposition of 
cos k x cos k y + (cosk x + cosk y ) representations with dif- 
ferent weights. As we move towards half filling, the hole 
pockets approach the Fermi level, the cos k x cos k y com- 
ponent is suppressed and the nodes of the gap move away 
from the electron pockets. 

Around half filling, 0.2 > /i > —0.3, the calculations 
converge to a SC gap of the form shown in Fig. [2]i which 
is of the nodeless s± type, cf. FigHb. In this regime, the 
FS consists of both electron and hole bands and the node 
lines lie almost halfway between them. This happens due 
to MD phenomena and the presence of interband cou- 
pling between the bands. The latter makes the SC gaps 
over the hole and electron FS coupled while the former 
ensures that maximum energy is gained by leaving the 
empty (no FS) space between these FS sheets ungapped. 
Such a mechanism results in a sign alternating s± solu- 
tion over the entire B.Z. The gap is nearly isotropic or 
may exhibit some arDOS driven momentum anisotropy 
(see Fig|3]) and corresponds to the nearly isotropic BCS- 
like gap reported in several experiments [El-Q • 

Upon further doping the system with holes, —0.35 > 
li > —0.55, the nodes of the s± gap intersect the outer 
hole FS (Fig 12b). Finally, at fi = -0.6 the FS and the 
arDOS distribution become as shown in Fig.[2f. The gap 
now is a superposition of s± and d x i_ y i symmetries with 
nodes at the low arDOS points of the FS. This rich struc- 
ture agrees with previous theoretical studies suggesting 
the near degeneracy of these two symmetries [17[ . Along 
with the pure c^-^-wave p-wave and s± solutions, this 
gap structure offers a possible explanation to NMR and 
penetration depth reports of nodal superconductivity in 
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FIG. 4. (Color online) Color grade plots of characteristic 
superconducting gap amplitudes in the folded BZ calculated 
self-consistently, when moving from the electron to the hole 
doped regime for q c = 7r/6. The panels a), b), c), and d) cor- 
respond to fi = 1.15,0.8,-0.1, and —0.6, respectively. The 
corresponding FS shapes are depicted in Fig. [2] Thin blue 
lines indicate the locations where the obtained SC gap van- 
ishes. In the color grade, red color illustrates positive and 
blue a negative sign of the maximal gap. 



both electron and hole doped compounds [H, [Toj| . No- 
tice that our theory predicts transitions from nodeless to 
nodal SC when going from the underdoped to the over- 
doped regime, (either with hole or electron doping). Ex- 
actly such a behavior has recently been found experimen- 
tally m. 

In conclusion, we have self-consistently computed the 
BCS gap of electron and hole doped iron pnictides. Our 
study offers a scenario for understanding the gap sym- 
metry of the superconducting state in these materials. 
We have shown that small-q electron-phonon scattering 
can accommodate conflicting experimental results pro- 
ducing a gap that changes sign while being compatible 
with the observed isotope effect. When both electron 
and hole FS sheets are present, we obtain within this 
purely phononic mechanism as self-consistent solution 
the fully gapped s± state, which was so far considered 
as evidence for a spin-fluctuations mechanism. Depend- 
ing on doping, the computed gap symmetries vary from 
nodeless to nodal in agreement with recent observations. 
Our solutions include s± + d-wave in the hole doped side 
and d-wave, as well as spin- triplet gap functions in the 
electron doped side. The present results establish that 
the small-q phonon picture for the pairing mechanism 
in the iron pnictides can explain the conflicting obser- 
vations and, hence, is the plausible one. Understanding 
the competition of SC with spin density waves as well as 
magnetoelastic effects when forward processes dominate 
the electron-phonon interaction could provide further in- 
sight. 
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